2005 International Linear Collider Workshop - Stanford, U.S.A. 



Two-loop contributions to electroweak precision observables in the MSSM 

J. Haestier, D. Stockinger, G. Weiglein 

IPPP, University of Durham, Durham DH1 3LE, U.K. 

S. Heinemeyer 

CERN, TH Division, Dept. of Physics, 1211 Geneva 23, Switzerland 



The evaluation of the two-loop MSSM-contributions to the electroweak precision observables M\y and sin 2 9 e g at 
O(o 2 ), 0(«t«(,), 0(a 2 ) is presented. These contributions enter via Ap, and it is explained in detail how one can 
retain the true, non-vanishing value of the MSSM Higgs boson mass in spite of using the gauge-less limit in the 
calculation. The numerical results can be sizeable, in particular for strong squark mixing. By comparing the results 
' in the on-shell and DR renormalization schemes, the remaining theoretical uncertainty is found to be small. 
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*-*\ 1. INTRODUCTION 

| Electroweak precision observables (EWPO) like Mw or the effective leptonic weak mixing angle sin 2 c e are inti- 
J> mately related to the quantum structure of the electroweak interactions. The experimental determination of these 
quantities, obtained at LEP and Tevatron, has a precision of better than one per-mille: SMw = 34 MeV (0.04%) 
and 6sm 2 6 e fi = 16 x 10~ 5 (0.07%) [1]. In the future, these accuracies will improve, and at the GigaZ option of a 
[ linear e + e~ collider a precision of 5M\y = 7 MeV [2, 3] and 5 sin 2 ^efi = 1-3 x 10~ 5 [3, 4] can be achieved. 

These precise measurements constitute tests of the quantum level of the Standard Model (SM) that are sensitive 
even to two-loop effects. The corresponding theoretical evaluation of the SM predictions up to the two-loop level is 
quite advanced, see in particular Ref. [5] for the most recent developments. On the other hand, the measurements of 
O the EWPO can be used to discriminate between different models of electroweak interactions and to derive constraints 
on unknown parameters. The comparison of the SM and its minimal supersymmetric extension, the MSSM, is 
particularly interesting since the MSSM agrees with all precision data at least as well as the SM, in some cases even 
better. 
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It is therefore highly desirable to know the MSSM predictions for the EWPO with a precision that matches the 
one of the SM and the experiments. The comparison of the SM and the MSSM predictions with the data could then 
lead to precise constraints e.g. on masses of supersymmetric particles. As a step towards this goal we have evaluated 
[6] the two-loop MSSM-corrections to the EWPO that enter via Ap at 0(a 2 ), 0(a t a b ), C(a 2 ), where 

^)_^(0) 

P Mf M w K ' 

in terms of the Z and W self energies Yiz.w- These are leading two- loop contributions involving the top and bottom 
Yukawa couplings and come from three classes of diagrams as shown in Fig. 1. These contributions to Ap induce 
universal two-loop corrections to the EWPO as follows (with 1 — s w = c w = M w jM\ ): 

„2 „2 



Mw C w A ^ \ c ,;„2 fl ^ C W S W 



SM W * 2 w ^ 2 Ap, 5sm'8 eS *- 2 w w 2 Ap. (2) 



The previously known two-loop contributions to EWPO in the MSSM comprise only QCD and SUSY-QCD corrections 
[7] and the C(a 2 ), 0(a t ab), C(a 2 ) corrections for the class (q) in Fig. 1 [8]. Obviously, the diagrams of class (q) 
contain no supersymmetric particles. The contributions from classes (q), (H) considered here in addition can be 
expected to be important for squark/Higgsino masses in the electroweak range and to have a more pronounced 
dependence on the MSSM parameters. A similar situation was found for the case of (g — 2) M , where the two-loop 
contributions from squark- Higgs diagrams are more important than the ones from quark- Higgs diagrams [9]. 
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Figure 1: Sample diagrams for the three classes of contributions to Ap considered here: (q) i/6-quark loop with Higgs exchange, 
(q) i/&-squark loop with Higgs exchange, (H) mixed quark/squark loop with Higgsino exchange. 

In the following we will first discuss in detail the renormalization and the restrictions imposed by approximating 
the EWPO by Ap as in cq. (2). It turns out that a strict calculation of Ap at 0(af), 0(a t (Xb), 0(a 2 ) would imply a 
vanishing MSSM Higgs boson mass, Mh — 0, which would lead to a bad approximation for the EWPO. We will show 
that it is possible to improve the approximation by taking into account the true value of Mh essentially everywhere. 
Finally we will discuss the numerical results and give an account of the remaining theoretical uncertainty based on 
the dependence on the renormalization scheme. 

2. THE ROLE OF THE GAUGE-LESS LIMIT AND THE LIGHT HIGGS BOSON MASS M h 
2.1. Gauge-less limit and M h = 

Ap as defined in eq. (1) constitutes a part of the loop corrections to the EWPO but is no observable. By itself, 
Ap is only UV-finite and gauge independent if corrections are considered that enter the EWPO only through vector 
boson self energies in the form (2) (as opposed to e.g. vertex or box diagrams). This is the case for the 0(a^), 
0(atab), O(o%) corrections. 

Taking into account strictly only terms of 0{pq ), 0(atab), 0(a^) corresponds in particular to neglecting the gauge 
coupling a and thus to taking the gauge-less limit a — > 0. In this limit, Mz,w while the ratio cw = Mw/Mz is 
fixed. 

In the SM the gauge- less limit is a reasonable approximation since a<at. In the MSSM, however, the gauge- less 
limit has side-effects in the Higgs sector since supersymmetry relates the Higgs self couplings to gauge couplings. At 
tree-level, the gauge-less limit implies the relations 

M h = 0, (3a) 
M H ± = M H = M A , sin a = — cos/3, cos a = sin/3 (3b) 

for the light and heavy CP-even, the CP-odd and the charged Higgs boson masses Mh,n, Ma, M h ± and the angles 
a, f3 where tan (3 — 1)2/1)1, the ratio of the two vacuum expectation values. 

The particularly troublesome of these relations is the masslessness of the light Higgs boson, Mh = 0. In the SM, 
where the Higgs boson mass is a free parameter even in the gauge-less limit, one knows that the result for Ap at 
O(o.f), 0(a t ceb), 0(al) is proportional to a factor [10, 11] 

19 - 2tt 2 for M H = 0, 19 - 2tt 2 + f(M H ) for M H ^ 0. (4) 

For typical values of Mh =O(100 GeV), this factor is about an order of magnitude larger than for Mh — 0, and the 
result for Mh = leads to a bad approximation for the EWPO. 

Due to the similar structure of the SM and MSSM diagrams it can be expected that taking the gauge-less limit 
relation Mh = literally would also lead to a bad approximation for the EWPO and should therefore be avoided. 
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Indeed, in Ref. [8] it was observed that the class (q) contributions to Ap in the MSSM are already UV-finite if all 
relations in eq. (3b) but not = are employed. In the following, we give an explanation of this result and show 
how it extends to the contributions of classes (q), (H). 



2.2. Comparison of the MSSM and a general two-Higgs doublet model 

In order to study these questions it is very useful to regard the MSSM as a special case of a more general two- 
Higgs-doublct model (2HDM) without supersymmetry relations for the couplings. Ap can be calculated at O(af), 
0(a t ab), O(al) in the gauge-less limit in both models, but in a general 2HDM, the gauge-less limit does not enforce 
any of the relations in eq. (3). 

Comparing first the contributions from class (q) in the MSSM and the 2HDM, we find that the corresponding 
two-loop diagrams and the counterterm contributions from the top/bottom sector are identical. The only difference 
concerns the Higgs sector counterterm contributions. In a general 2HDM, they can be derived from the one-loop 
expression (F is a symmetric function satisfying F (x,x) — 0, d x F (x,y)\ y=x = 0) 

Ap^hdm « [Fo(M 2 h± , Mi„) + 8 }_ a (F (M 2 H± , Ml) - F (M 2 A „, M%)) + c 2 _ a (F (M 2 H± , M 2 ) - F (M 2 ,M 2 )) ] (5) 

by performing the renormalization transformation Mh — > Mh + 5Mh etc. In the MSSM, eq. (5) and the Higgs sector 
counterterms vanish because of the gauge- less limit relations (3). Thus the class (q) contributions can be decomposed 

Ap^ Loop + Ap[%_ cta + Ap<£ cts in the 2HDM and as Ap^ Loop + Ap[%_ cta in the MSSM. The 2HDM result is 
UV-finite for all values of the Higgs sector parameters. 

From this comparison we find that the MSSM result is not only UV-finite if (3) is used but more generally if 
Apn-cts as derived from (5) is finite. From eq. (5) we can explicitly read off that Ap^_ cts = already if only the 
relations (3b) are used. This explains the observation made in Rcf. [8]. The MSSM result Ap^Loop + Ap[^_ cts with 
the relations (3b) corresponds to a certain special case of the general 2HDM calculation. Therefore it is finite even 
for the true, non-vanishing value of Mh- 

For the class (q, H) contributions we obtain similar decompositions 

P2HDM — ^2-Loop + ^t/b-cts ' ^t/b-cts + ^H-cts' W 
-Loop + ^Pt/b-cts + ^t/b-cts 



^PMSSM — ^^-Loop + ^Pt/b-cts + ^P* /h-rt« ■ {') 



In this case, even if the relations (3b) arc used such that Ap^l^ s vanishes, there is still a difference between the 
2HDM and the MSSM result because the t/b sector counterterms differ as indicated by the superscripts OS4 and 
OS3. 

In the MSSM, supersymmetry in conjunction with SU(2) gauge invariance correlates the four sfermion masses 
m ti 2 ' m b 1 2 ' Therefore only three of them can be renormalized independently. We choose to renormalize 2 , and 
mg 2 independently by on-shell conditions (OS3). Then the fourth renormalization constant Sm^ is determined as a 
function of the other three in the MSSM, while in the 2HDM all four sfermion masses can be defined independently 
by on-shell conditions (OS4). In the simple case of vanishing mixing between left- and right-handed sfcrmions, 
m bi = m b L an< ^ wc obtain the following different results for 6m? : 

8mlJos3 = 8m 2 L +5m 2 b -5m 2 , 5m?Jos4 = ^(mj-J. (8) 

The difference between the sfermion sector counterterms in the 2HDM and the MSSM is thus contained in the mass 
shift 

Am l = H L ) - ( Sm l + Sm b - Sm t) ■ (9) 

It turns out that this mass shift is only UV-finite if all gauge-less limit relations are used, including Mh = 0. 
Accordingly, the MSSM result (7) is only finite if all relations in (3) are employed. The 2HDM result with on-shell 
renormalization of all four sfermion masses is of course UV-finite for all choices of Mh ■ 
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Figure 2: Ap^ and Ap (,3 " 2HDM) as functions of M h for M SUSY = 400 GeV, A t = 1250 GeV, tan/3 = 2, fj, = 500 GeV and 
M A =150 GeV (short-dashed), 300 GeV (full line), 1000 GeV (long-dashed). This scenario involves a light stop with m t ~ = 120 
GeV. 

2.3. Incorporating the Higgs boson mass into the MSSM result 

The best way to take into account the non- vanishing value of M h as much as possible is to consider the combination 

^il(M h = 0) + (Ap%£l(M h ) Apg*i(M h = 0)) , (10) 

where the relations (3b) are used everywhere. As explained above, all three terms are individually UV-finite. Since 
the difference between Ap^ssM and A^hdm ^ s connnc d to the t/b counterterms, the combination (10) can be written 
as 

A/4-Xop(M.) + ±p?£Lm) + A^™ 4) (M.) + [Apg"£° S3) (M, = 0) - A P ^\M h = 0)] . (11) 

The first three terms correspond to the MSSM calculation where Mh is set to its true, non-vanishing value but 
where all sfermion masses are renormalized by on-shell conditions instead of using the mass relation imposed by 
supersymmetry. The term in the square brackets is proportional to the mass shift Am? (Mh = 0) that restores the 
necessary supersymmetry mass relation. It is only here that Mh = has to be employed. 



3. NUMERICAL RESULTS 

3.1. Results for different values of the supersymmetry parameters 

Fig. 2 demonstrates that it is in general important to take into account the true value of the Higgs boson mass 
Mh- In the left and right panels of Fig. 2, Ap^ and Ap^ BM , which is the only M^-dependent term in cq. (10), are 
shown as functions of Mh in a scenario with a light stop (see caption). For both the fermion and the sfermion loop 
contributions the difference of setting Mh — or Mh =O(100 GeV) amounts to more than 10~ 4 . 

In the remainder we focus on the numerical results from the classes (q, H) (for class (g) see Ref. [8]). In Figs. 3, 4 
the results for Ap as defined in eqs. (10,11) are shown, split up into the contributions for the individual classes. We 
choose several representative values for the supersymmetry parameters as described in the captions, basically always 
starting from the SPSla scenario [12] and varying one or two of the parameters. In the class (q) contributions the 
Higgs boson mass Mh is set to either 120 GeV, which is a good approximation for the true, loop-corrected value of 
Mh, or to zero; the class (H) contributions are M^-independent. 

According to eq. (2) a contribution of Ap = 10~ 4 leads to shifts 

SM W = 6 MeV, 5 sin 2 6» cff = -3 x 10~ 5 . (12) 

We find that contributions of this order of magnitude are possible, in particular for strong sfermion mixing (large 
values of ji or A t in Fig. 4) but also for a light common sfermion mass parameter MgusY as in Fig. 3. 
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Figure 3: Results for Ap W) (M ft = 120 GeV) (l,brown), Ap (H) (2,green), and, for comparison, Ap^(M h = 0) (3,red), shown 
as functions of the common sfermion mass -Msusy. In the left panel, the A t ,t = 0; in the right panel, A tt b are chosen such 
that the ratios Msusy : At : Ab are as in the SPSla scenario. The remaining supersymmetry parameters are always set to the 
values of the SPSla scenario. 
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Figure 4: Results for Ap {i) (M h = 120 GeV) (l,brown), Ap (H) (2,green), and, for comparison, Ap {i) (M h = 0) (3,red). In the 
left panel, the results are shown as functions of the /i-parameter for tan /3 = 40, such that for large p one sbottom becomes 
light; in the right panel, the results are shown as functions of At- The remaining supersymmetry parameters are always set to 
the values of the SPSla scenario. 



3.2. Different renormalization schemes and estimate of remaining theoretical uncertainty 

So far we have chosen the on-shell renormalization scheme for the independent sfermion masses 2 and . 
In Fig. 5 the on-shell results are compared with the results in the DR scheme, where the counterterms for the soft 
supersymmetry breaking parameters in the sfermion mass matrices are defined as pure divergences. This DR scheme 
implies non-vanishing finite parts of the sfermion mass counterterms, which read, in the case of vanishing left-right 
mixing, Sm 2 f |fi n _ pa rt — 5m 2 f |fl n _ pa rt- The comparison of the two renormalization schemes is interesting in order to 
assess the numerical stability of the result and the intrinsic theoretical uncertainty of the two-loop contributions. 

We find that the difference between the on-shell and DR results for the sfermion loop contributions at the one- 
loop level is of the order 10~ 5 to 10 -4 for large sfermion mixing. At the two-loop level this renormalization-scheme 
dependence is significantly reduced to well below 10 -5 . We have checked that this result is general and not restricted 
to the particular parameter choice in Fig. 5. 
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Figure 5: The sfermion loop one-loop result Ap' 1L ^ (dashed) and the sum Ap (1L) + Ap^ H) (full lines) in the on-shell (blue) 
and the DR (red) renormalization scheme. The parameters are as in the right panel of Fig. 4. The left panel shows the full 
results, where the differences between the curves are hardly visible; the right panel shows the differences of the results to 
a P (1L) Idr. i- e - C 1 ) shows A/9 (1L) | os - Apgy, (2,3) show (Ap (1L) + Ap {9 ' 6) }aK,os ~ A P^m ■ For simplicity, M h = is used 
here, but the conclusions do not change for non- vanishing Mh- 

In conclusion, we have evaluated the 0(a|), 0(a t ah) : C(a 2 ) corrections to Ap and thus to the EWPO Mw and 
sin 2 # c ff in the MSSM. Although the gauge-less limit is necessary and leads to the tree-level relation Mh = 0, we have 
shown that the true MSSM Higgs boson mass can be taken into account. The numerical values of the class (q, H) 
contributions (10), (11) to Ap can amount to 10 -4 , corresponding to shifts 5Mw = 6 MeV,5sin 2 8 c g — —3 x 10 -5 . 
The comparison of the two renormalization schemes shows that the inclusion of the two-loop result in the MSSM- 
prediction for Ap and the EWPO leads to a significantly improved accuracy. The residual theoretical uncertainty 
due to unknown three-loop corrections of 0{a^ b ) is well below the foreseen experimental resolution achievable at the 
GigaZ option of a linear e + e~ collider. 
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